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We review various combinatorial applications of field theoretical and matrix model 
approaches to equilibrium statistical physics involving the enumeration of fixed and random 
lattice model configurations. We show how the structures of the underlying lattices, in 
particular their color ability properties, become relevant when we consider hard-particles 
or fully-packed loop models on them. We show how a careful back-and-forth application 
of results of two-dimensional quantum gravity and matrix models allows to predict critical 
universality classes and consequently exact asymptotics for various numbers, counting in 
particular hard object configurations on fixed or random lattices and meanders. 
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1. Introduction 

Physics provides a source of alternative approaches to well-posed mathematical prob- 
lems, allowing sometimes for less rigorous but more predictive results, such as those dis- 
cussed in these notes. The reason is that the tools employed are often somewhat ill-defined 
from a purely mathematical point of view, although they are often validated by experi- 
ment. In every physicist's toolbox, field theory is certainly the most powerful one, as it 
basically allows for an effective description of the essential behavior and properties of the 
situations it models. 

Our main concern here will be with solving well-posed combinatorial problems by 
means of physical approaches and reasoning. This type of approach must have a sort of 
exotic flavor when observed from a purely mathematical point of view. The very principles 
it relies on borrow from the developments of field theory, such as renormalization group 
techniques applied to critical phenomena. Instead of trying to be completely exhaustive, 
which would go far beyond the scope of the present notes, we will rather try to extract 
the essentials and main lessons of the field-theoretical approaches. We will also present 
some alternative techniques such as matrix integrals allowing to address some of these 
combinatorial problems. 

The questions we will be considering here deal with fixed or random statistical two- 
dimensional lattice models, used to describe two-dimensional interfaces as well as discrete 
two-dimensional quantum gravity (2DQG). Combinatorially, these are simply problems of 
enumeration of configurations involving either decorated lattices or decorated graphs. More 
precisely, we will discuss here two classes of models for which the details of the underlying 
lattice are important: these models "feel" strongly their underlying discretized space. 
They are the hard-particle and fully-packed loop models. We will be mainly interested in 
understanding the critical behavior of these systems, namely in locating phase transitions 
in the space of parameters where thermodynamic quantities, characteristic of the collective 
behavior of the system at large size, become singular. As we will explain, these singularities 
fall into "universality classes" described by various field theories, and characterized by 
critical exponents governing their algebraic singularities. 

The first theory describes objects (particles) occupying the vertices of the lattice, with 
a nearest neighbor exclusion constraint that no two vertices connected by an edge may be 
simultaneously occupied. Despite the local character of this constraint, it allows for the 
particles to "feel" the details of the lattice. When for instance we increase the density 
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of occupation up to its maximum, particles will try to occupy every other vertex in the 
lattice, thus will feel for instance whether it is bipartite or not. 

The second theory describes loops occupying edges of the lattice, with the additional 
constraint that every vertex is visited by a loop. Again, the bipartite (or not) character of 
the lattice will prove crucial in the discussion of the critical behavior of the system. 

While solving the first theory will eventually lead to the exact enumeration of ran- 
dom planar graphs with hard-particles on them, solving the second one will lead to an 
asymptotic solution of the so-called "meander" problem of enumeration of the topologi- 
cally inequivalent configurations of a closed road crossing a river through a given number 
of bridges. The materials collected in this review article include in particular Refs. [Q and 
0, where random lattice hard-particle and meander problems were first solved (see also 
the review article 0). 

The paper is organized as follows. In Sect. 2, we collect a few definitions and results 
regarding the thermodynamics of two-dimensional statistical lattice models on fixed and 
random lattices, including in particular field-theoretical ones. Sect. 3 is devoted to the 
study of hard particles on various random lattices, in the form of planar graphs with fixed 
valence. We show in particular how another favorite physicist's tool, the matrix integral, 
may be used to completely solve various enumeration problems involving decorated graphs. 
When applied to hard particles, these allow for a complete understanding of the crystalliza- 
tion transition, namely a critical value of the density of occupation below which the system 
is fluid (particles occupy vertices quasi-randomly) and beyond which a crystal of particles 
forms. We show however that such a transition only exists on (fixed or random) lattices 
with a colorability property allowing for the existence of maximally occupied sublattices. 
The critical behavior is then traced back directly to the existence and structure of these 
sublattices. We show in particular that for bipartite lattices, the critical transition lies 
in the universality class of the critical 2D Ising model, thus proving on a random lattice 
some old fixed lattice conjectures concerning hard squares or hard triangles []5] |J. In 
Sect. 4, we address fully-packed loop models on fixed and random lattices. After recalling 
some known results in the fixed lattice case, we move to random lattices. Remarkably, just 
like in the hard-particle case, we find that according to whether the random lattices are 
bipartite or not, the critical universality class predicted is different. By suitably following 
a two-flavor fully-packed loop model from its fixed lattice version to its random lattice 
one, we arrive at the identification of the universality class for meanders, which allows 
in particular to predict the meander configuration exponent a. = 29 + , governing the 
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large n asymptotics of the meander number with 2n bridges M 2n ~ g~ 2n /n a . We gather 
a few concluding remarks in Sect. 5, as well as another simple prediction, as a challenge 
for mathematicians to prove. 

2. 2D Quantum Gravity and Asymptotic Graph Combinatorics 

2.1. 2D Lattice Models 

The archetypical systems studied in 2D equilibrium statistical mechanics are lattice 
models, whose configurations are defined as maps a from, say the vertices i of a 2D lattice 
£ (or more precisely a finite connected subset thereof say of rectangular shape of size LxT, 
denoted by £l,t, and with definite, say periodic boundary conditions) to a target space 
T usually real. The model is further defined through an energy functional E(a), and each 
configuration is attached a statistical weight proportional to e~^ E ^ a \ where f3 = 1/(/cbT), 
ks the Boltzmann constant and T the temperature. One is usually interested in the 
properties of global functions such as the partition function 

where the sum extends over the set of maps from Cl,t to T. This is precisely the normal- 
ization of the probability p(a) = e~ f3E ( a > /Z LiT of the configuration a. This probability 
weight allows to answer various questions regarding correlations in the system, through 
the corresponding expectation value denoted generically by (...) = J2 a ...p(a). 

As a simple illustration, the Ising model has T = {— 1, +1}, and the images a(i) are 
used for instance to describe the spins (magnetic moments) in metals. The functional 
energy usually incorporates some information on interactions within the system or with 
some external fields or forces. For the Ising model, it reads E{a) = —JJ2(ij) a (i) a (j) ~ 
H ^2 1 <r(i), where the first sum extends over nearest neighboring pairs of vertices (i, j) and 
expresses ferro- (J > 0) or antiferro- (J < 0) magnetic interactions, while the second 
sum expresses the coupling of the system to an external magnetic field H. In general, the 
energy functional depends on external parameters (such as (3 J and j3H in the Ising case). 

The archetypical questions one tries to answer involve the thermodynamic limit in 
which the system becomes large, say L, T — > oo with L/T finite, and on its dependence on 
the external parameters. In this limit, thermodynamic functions such as Z = Hu^Z^t) 1 ^ 
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may develop singularities when the parameters approach some critical values, correspond- 
ing to phase transitions and critical phenomena. Of particular interest are those in which 
a divergence occurs for the length characteristic of the effective range of interactions, the 
correlation length. The corresponding critical singularities fall into so-called universality 
classes, characterized by scaling exponents which govern the leading algebraic behavior 
of the singular part of the thermodynamic functions of the system. For instance, in the 
ferromagnetic Ising model with vanishing magnetic field, one defines a "thermal" expo- 
nent a, characteristic of the leading singularity of the free energy F = LogZ when the 
temperature T approaches a critical value T c : 

F| sing oc \T-T c \ 2 ~ a (2.2) 

Similarly, at T = T c , we also get a singular behavior for F when we switch on a small 
magnetic field H, with 

F| sing oc {HI 2 -? (2.3) 

which defines the "magnetic susceptibility" exponent (3. The exponents a, (3 are expected 
to be more universal than the actual value of the critical temperature T c : indeed, one finds 
for instance that T c depends on the precise type of lattice on which the model is defined, 
while a, (3 only depend on its dimensionality. 

In two dimensions, it has been shown that under some assumption of locality of the 
interactions, critical singularities corresponding to the divergence of the correlation length 
fall into universality classes described by two-dimensional conformal field theories (CFT), 
upon taking the continuum critical limit of the models, in which we let at the same time the 
size of the system tend to infinity and the lattice spacing tend to zero while parameters tend 
to critical values, so that the lattice becomes the two-dimensional plane (or a higher genus 
Riemann surface according to boundary conditions) , and configuration maps become fields. 
Conformal invariance of the resulting field theory is directly linked to the divergence of the 
correlation length which implies in particular local scale invariance of the theory. Two- 
dimensional CFT's have been extensively studied and partly classified, providing us with 
a sort of Mendeleiev table of universality classes of critical phenomena in two dimensions, 
in the form of a list of critical exponents (see the text and references therein). More 
precisely, CFT's are characterized by their central charge c m = 1 — 6/(m(m + 1)) and 
the conformal dimensions of their operators h r>s = (((m + l)r — ms) 2 — l)/(4m(m + 1)) 
respectively the central extension and highest weights of the corresponding representations 
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of the Virasoro algebra, m some complex parameter, and r, s some positive integers, further 
restricted to a finite set in the case of minimal models, where m is rational. The conformal 
dimensions govern the fall-off with distance of correlation functions. For instance we have 
for a conformal operator $^ of dimension h: z)&h(0, 0)) = l/|z| 4/l . Upon identifying 

the operators of the theory as generating perturbations of the various parameters away 
from their critical values, one may relate the corresponding critical exponents to their 
conformal dimension through! 

1 - 2h . AS 

«=— (2.4) 

For instance, the critical 2D Ising model corresponds to a CFT with m = 3, i.e. c = 1/2 
and the thermal operator associated with perturbations in the temperature has = 1/2, 
so that the thermal exponent is a = (in this case, the singularity is weaker and involves 
the logarithm of \T—T C \). Besides, the magnetic "spin" operator, the continuous version of 
the map ex, has conformal dimension h,2,2 = 1/16, which leads to the magnetic susceptibility 
exponent (3 = 14/15. 

A number of classifications of CFT's are now available, the most remarkable of which is 
probably that of "minimal models" with central charges c < 1, namely CFT's with a finite 
number of conformal primary fields ||. They turn out to have central charges c = c(p, q) 
1 — 6(p — q) 2 /(pq), indexed by two coprime positive integers p and q. The constraint of 
unitarity (positivity of correlators, and in particular positivity of the probabilities p(cr)) 
further restricts q and p to obey \p — q\ = 1. In this language, the Ising model corresponds 
to p = 3, q = 4. 

2.2. 2D Random Lattice Models 

In this work, we consider random lattice models, defined in precisely the same way as 
ordinary lattice models, but by replacing the underlying lattice by a somewhat arbitrary 
possibly disconnected tessellation of the plane (or higher genus surface). The tessellation 
becomes therefore part of the configuration to be summed over (see Fig. 1 for an illustra- 
tion). This was introduced as a discrete model of 2D quantum gravity (2DQG), namely to 
describe the coupling of some statistical "matter" models to the quantum fluctuations of 



2 This is done by dimensional analysis. Writing the corresponding perturbation of the con- 
formal action functional as (t — t c ) J d 2 z&h(z, z), and performing a local rescaling z — > Xz, 
we find that (t — t c ) ~ \ 2h - 2 as ^ e ac tion is dimensionless. The thermodynamic free energy 
F = lim -^Log Zl,t however is per unit of area, hence (t — t c ) 2 ~ a ~ A~ 2 , hence a = 1/(1 — h). 
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Fig. 1: Schematic representation of how we go from a 2D statistical lattice 
model (defined here on a finite portion of the triangular lattice) to a random 
lattice model, in which the lattice is replaced by arbitrary tessellations of 
space (here a random triangulation tessellating a genus one surface) . 

the underlying "space", represented by random lattices or tessellations (see e.g. Ref. @ 
for a review and more references). For each tessellation O, we associate a statistical weight 
directly borrowed from the Einstein action in 2D, namely a weight N x (®'g A (®> where x 
and A are respectively the Euler characteristic and area of the tessellation (respectively 
measured via x = #faces— #edges+#vertices and A = Utiles), and where iV and g are 
the discrete counterparts of the Newton constant and cosmo logical constant. For each 
tessellation O together with a spin configuration a : — > T on it, we also have a weight 
e -pE{a) anc j we mus t divide by the order of the automorphism group \Aut(@, a)\. This 
leads for instance to the discretized partition function of any 2D statistical lattice model 
coupled to 2D quantum gravity 

Z= JV* ( «V (e) E e- m<T) /\Aut(e,a)\ (2.5) 

teSSel. config. 

cr : @ — y T 

Note that the parameter N allows to isolate the contributions of tessellations of fixed genus. 
In the following we will be mainly interested in the genus zero contributions, obtained by 
letting iV — > oo. Note also that the free energy F = LogZ selects only the connected 
tessellations in the sum ( |2.5| ). 

Like in the fixed lattice case, we are interested in the thermodynamic limit of the 
system, in which say the average area or some related quantity diverges, ensuring that 
the dominant contributions to Z e.g. come from large tessellations. This is guaranteed 
in general by the existence of a critical value g c of the cosmological constant g at which 
such singularities take place. This value is a priori a function of the type of random 
lattices we sum over as well as of the various matter parameters. We may now attain 
interesting critical points by also letting the matter parameters approach critical values, 
a priori distinct from those on fixed lattices. The result is well described by the coupling 
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of the corresponding CFT's with quantum fluctuations of space, namely by letting the 
metric of the underlying 2D space fluctuate. Such fluctuations may be represented in the 
conformal gauge by yet another field theory, the Liouville field theory, which is coupled to 
the matter CFT. This field theoretical setting allows for a complete understanding of the 
various critical exponents occurring at these critical points ||21|| . For instance, one defines 
the (genus zero) string susceptibility exponent 7 s t r as the exponent a associated to the 
cosmological constant singularity, namely by writing the singularity of the free energy as 

rising ^(g c -9) 2 - l3tr (2.6) 

In the case of coupling of a matter theory with central charge c to 2DQG, one has the 
exact relation 

. Ur = C -W(l- C )(25-c) p 7) 

In the case of "pure gravity", namely when the matter is trivial and has c = 0, we get 
Istr — — l/2, while for the critical Ising model with c = 1/2 we have 7 s t r = —1/3. 

Upon coupling to gravity the operators of the CFT (&h(z, z)) get "dressed" by gravity 
($/j, — » = ^a) and acquire dressed dimensions A, given similarly by 

A = Vl-c + 24 ft -VT^ 

V25 - c - ^/T^c K ' 

As opposed to the fixed lattice case, where conformal dimensions govern the fall-off of 
correlation functions of operators with distance, the dressed operators of quantum gravity 
do not feel distances, as their position is integrated over the surfaces, but rather only feel 
changes of area at fixed genus. More precisely the general genus zero correlators behave 
in the vicinity of g c pi] as 



These results may be easily translated into the large (but fixed) area A behavior of 
the various thermodynamic quantities, upon performing a Laplace transform, which selects 
the coefficient of g A in the various expansions. Let Fa denote the partition function for 
connected tessellations of genus zero and area A, we have 

q~ A 

F A ~ ~ir (2-10) 
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while if (...)a denotes any genus zero correlator at fixed area, we have 

(*A 1 *A 2 ...*A„>A ~ — ^ (2.11) 

In [21], all these formulas were also generalized to higher genus as well. 

As an illustration of eqn. ( 2.10 ), recall that the number of quadrangulations of the 
sphere with A square tiles and with a marked edge (dually equal to the number of rooted 
tetravalent planar maps with A vertices) reads |n 



3 A /2A\ 12 A 

Na= 2(A + 1)(A + 2){a)~AW (2 ' 12) 

Noting that the rooting simply amounts to Na oc AFa, the asymptotics Q2.12Q correspond 
to g c = 1/12 and 7 s ^ r = —1/2, hence c = 0. This is one of the various ways to attain 
the universality class of pure gravity, namely by summing over bare tessellations without 
matter on them. 

So far so good, basically all asymptotic combinatorics problems involving planar (or 
more generally fixed genus) graphs seem to be solved via eqns. (|2.10 ) and ( |2.11|) or their 
higher genus generalizations, provided one is able to identify the central charge c of the 
underlying CFT. This latter step however may prove to be quite involved. In fact, the 
aim of this note is to clarify a large class of cases for which the naive application of these 
formulas would lead to the wrong result. These are statistical models whose definition 
strongly relies on the structure of the underlying (fixed or random) lattice. This is the 
case for hard particles, as well as for fully-packed loops, the two main subjects of this note. 



3. Hard Particles 

The first class of problems we will be discussing is that of hard particles occupying or 
not the vertices of the fixed or random lattices, and subject to the "hardness" or particle- 
exclusion constraint that no two adjacent vertices may be simultaneously occupied. One 
usually attaches an activity z per occupied vertex. In this section, we first recall the 
fixed lattice results, and display the exact or predicted critical behavior of the model. 
Generically it undergoes a crystallization transition between a disordered phase of low- 
density of occupation and a crystalline phase of maximal occupancy. We next turn to 
the exact solution of the same model on some random lattice, only to discover that the 
crystalline transition is wiped out by the sum over lattices. To solve this puzzle, we simply 
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notice that the transition can take place only if the random lattices allow for the existence 
of some generic crystalline configurations. We finally show that the model, when defined 
on random bipartite lattices, has a crystalline transition in the universality class of the 
two-dimensional critical Ising model coupled to 2DQG. 



3.1. Fixed lattice results and conjectures 
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Fig. 2: Typical configurations of hard particles on the (a) hexagonal, 
(b) square and (c) triangular lattices. We have also represented in grey 
the excluded area translating the nearest neighbor exclusion into an non- 
overlapping constraint of tiles. Their respective shapes are triangles, squares 
and hexagons, hence the names hard triangles, squares or hexagons. 



Hard particle models have been extensively studied in the physics literature, under 
various names: hard-core lattice gas, nearest-neighbor exclusion models, etc ... [0 @ 



||iq| |TT |. Among the host of results and conjectures, let us mention the three cases of the 
hexagonal, square and triangular lattices. The corresponding models are also often referred 
to as hard triangle, square or hexagon models, as in the dual picture an occupied vertex 
may be replaced by some tile whose shape transforms the hard-particle constraint into a no- 
overlap condition between tiles (see Fig. 2 for an illustration in the hexagonal (a), square 
(b) and triangular (c) cases, respectively corresponding to hard triangles, squares and 
hexagons). In the three cases, one expects in the thermodynamic limit a crystallization 
transition to take place at some critical value z = z+ > between a low-z disordered 
phase and a high-z crystalline phase in which, when z — > oo, only one of a few available 
"groundstates" are realized, namely some maximally occupied configurations, in which the 
particles occupy a sublattice of the original one. For the hexagonal and square cases, where 
the lattice is bipartite, there are exactly two such crystalline groundstates, as opposed to 
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the triangular case, where there are three. In addition, one expects as well some "non- 
physical" critical point at some other value z = Z- < 0, independent of the precise details 
of the lattice, and generically described by the Yang-Lee edge singularity |L2| (of the Ising 
model in the presence of a critical imaginary electric field). 

The triangular lattice case is particularly interesting, as it proved to be exactly solv- 
able by Baxter [[TO], who explicitly worked out various critical exponents characterizing 



the crystallization transition. He found the two critical points z± = ( 1± 2 V ^' ) 5 ' w hh re- 
spective thermal exponents ct_ = 7/6 and a+ = 1/3 governing the singularities of the 
thermodynamic free energy f\ s %ng ~ \ z ~ Zi\ 2 ~ ai , % = ±. These correspond respectively 
to the minimal CFT's with (p, q) = (2, 5) and (p, q) = (5, 6) respectively, namely with 
central charge c(2, 5) = —22/5 and c(5,6) = 4/5. These are the universality classes of 
the 2D Yang-lee edge singularity |13[ and of the 2D critical three-state Potts model (a 
generalization of the Ising model with three spins instead of two). 

As to the two other cases of hard squares and triangles, no exact solution is known 
to this day. The powerful Corner Transfer Matrix method of Baxter has allowed for many 
exact series expansions for thermodynamic quantities up to quite large orders [|10| [p 



eventually leading to the very reasonable conjecture that in these two bipartite cases the 
crystallization transition point lies in the universality class of the critical Ising model, at 
some critical value z+ depending on the lattice, while the other critical point still lies in 
the Yang-Lee edge singularity class, at some z_ also depending on the lattice. 

Beyond numerical evidence accumulated so far, we will propose here a "gravitational" 
proof of the identification of these universality classes via their random lattice version. 

3.2. Matrix models as combinatorial tools 

Matrix models are remarkably efficient tools to realize the summation over (possibly 
decorated) tessellations of given genus of the type of ( |2.5|) (see e.g. [[L4| and |15j for re- 



views). The basic idea is to replace the task of performing say a Gaussian integration 
over an iV x N Hermitian matrix M by the drawing of pictures, called "Feynmann dia- 
grams" , weighted accordingly in such a way that the weighted sum over diagrams equals 
the Gaussian integral. 

As an archetypical example of such a thing, let us display the quartic one-matrix 
model, with partition function 

Z(g) = J c!Me- NT ^-^ (3.1) 
10 
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M u Tr(M 4 ) <M U M tl > 

Fig. 3: The pictorial representation for matrix elements Mij (double half- 
edges carrying matrix indices), terms Tr(M 4 ) (a collection of four double 
half-edges connected at a vertex, with matrix indices conserved along the 
oriented lines), and the propagators used to glue pairs of double half-edges 
into edges in the fatgraphs representing the matrix integral. 

where dM stands for the standard Haar measure, normalized in such a way that Z(g = 
0) = 1. Z(g) is to be understood as a formal power series of g, the coefficients of which 
must be computed by performing Gaussian matrix integrals. By virtue of Wick's theorem, 
this is readily done (for the coefficient Z n of g n ) by computing 

J nl 4 

^ £ n 

pairs (ij),(kl) ij,kl 

where the sum extends over all decompositions of the matrix elements into pairs and the 
notation (M y M fcI ) = jjO~u&j,k stands for the "propagators" / dMe~ NTr ( M2 M i:j M kl . In 
turn, the value of the propagator allows to devise a pictorial representation for the sum 
in ( p.2| ). Representing a matrix element My by a double half-edge as in Fig. 3, where 
each oriented line carries a matrix index running from 1 to N, we may express Tr(M ) as 
in Fig. 3, and interpret the sum in ( |3.2|) as that over all possibly disconnected fatgraphs 
obtained by gluing half-edges into edges (propagators) . Each such graph receives a weight 
1/N per edge (from the propagators), Ng per vertex, and N per remaining oriented loop 
carrying an index summed from 1 to N. This gives an overall weight g A j\f F ~ E + A for each 
possibly disconnected graph (A = ^vertices, F = #faces and E = hedges), as well as 
an overall rational factor expressed as l/(4 n n!), which is nothing but the inverse of the 
order of the automorphism group of the corresponding fatgraph. Finally upon taking the 
logarithm of Z(g), one ends up with the desired sum over connected tessellations © (made 
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(3.2) 



of squares here, dual to the tetravalent vertices of the matrix model), and with a weight 
]yx(&) g A ( e ) f\Aut(Q)\. This particular example therefore corresponds to pure gravity, 
in the form of random quadrangulations. Note that the size N of the matrix serves as 
Newton constant in this approach, while the cosmological constant g is a feature of the 
potential (term in the non-gaussian part of the exponential weight). 

More generally, we may engineer matrix models to suit our needs, either by considering 
a more general potential allowing for weighted vertices of arbitrary valencies, or/and by 
considering multi-matrix integrals of the form 

Z = f dA^dA*- 1 ™^**-^ v ^ (3.3) 

where A denotes the column vector of matrices (A±, ...,AkY and Q is a k x k quadratic 
form (with A t QA = Xa< a b<k QabA a Ab), and the Vi some polynomial potentials. Note 
that again we normalize the measure in such a way that Z = 1 if the Vi vanish. The 
diagrammatic interpretation of Z follows from the multi-matrix generalization of Wick's 
theorem expressing the multi-Gaussian matrix integral as a sum over pair decompositions 
of matrix elements of the integrand evaluated by using propagators, namely 

(f({(A a )i j })) Q = £ HdAaUAbhih (3.4) 

pairings pairs 

where (...)q = / Yl a dA a ...e~ NTr ^ At( ^ A ^ stands for the multi-Gaussian average. In par- 
ticular, propagators are easily computed in terms of the inverse of the quadratic form 
Q: 

((A a ) ij (A b ) kl ) Q = 6 -^(Q-\, b (3-5) 

The diagrammatic interpretation of Z is now clear: attach different "colors" a = 1, 2, k 
to the matrices. The potentials Vi(x) = ^ gi,mX m /m contain the list of allowed vertices 
of color i and valence m with their respective weights <7i, m . The quadratic form Q encodes 
via its inverse the allowed propagators, namely the list of possible edges linking these 
vertices to one another. The result is a sum over vertex-colored graphs with specific gluing 
rules. Note moreover that irrespectively of its color each vertex receives a weight N, each 
edge a weight 1/N and each loop of oriented index lines a factor N, which still produces 
the topological factor N x in the expansion of LogZ, x the Euler characteristic of the 
(decorated) connected graph. This allows in particular to extract the contribution from 
planar graphs by taking iV — > oo. 

In the next two sections, we use some of these models to solve the hard-particle 
problem on various random lattices. 
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Fig. 4: A typical configuration of hard particles on a tetravalent planar 
graph. The particles are represented by filled circles (•). They obey the 
nearest neighbor exclusion constraint that no two particles may be adjacent 
to the same edge. 

3.3. Hard objects on random tetravalent graphs: matrix model solution 

We now turn to the enumeration of hard-particle configurations on random tetravalent 
graphs, the natural fluctuating versions of the square lattice. For illustration we display 
in Fig. 4 a typical hard-particle configuration on a connected planar tetravalent graph. 
The occupied vertices (represented by filled circles) obey the nearest-neighbor exclusion 
(hardness) constraint that no two particles may be adjacent to the same edge. As explained 
before, we attach a weight z per particle and g per vertex. 

In view of previous section, we are now ready to write a two-matrix model for enu- 
merating such configurations, with partition function 

Z(g, z) = j dAdBe~ Ni: < AB -^-° zJ £-^ (3.6) 

where the measure is now normalized in such a way that Z(g = 0, z) = 1. The graphical 
interpretation (see Fig. 5) is a consequence of the two-matrix form of the Wick theorem, 
in which propagators are calculated using the Gaussian weig ht e -NTr(AB-A*/2)^ namely 
with (AijBki) = diidjk/N and (BijBki) = SuSjk/N while (AijAki) = as a particular case 
of eqn. ( |3.5| ). Representing the terms Tr(S 4 ) as empty tetravalent vertices and Tr(A ) as 
occupied ones, we see that the gluing of half-edges into edges (through the propagators) 
now incorporates the hard-particle exclusion rule that no two occupied vertices may be 
adjacent (vanishing of the (AA) propagator). Again, Log Z(g, z) is the generating function 
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Fig. 5: The pictorial representation for the matrix elements of the two- 
matrix model describing hard-particles on tetravalent fatgraphs. The occu- 
pied (resp. empty) vertices correspond to the confluence of four A (resp. B) 
matrix elements, represented as solid (resp. dashed) double half-lines. We 
have also represented the only two non-vanishing propagators allowing for 
gluing double half-edges into edges of the final graph. These obey the par- 
ticle exclusion rule that no two adjacent vertices of the final graph may be 
simultaneously occupied. This information is encoded in the quadratic piece 
of the potential of the two-matrix model. 

for connected fatgraphs with hard particles on them, and with a weight g per vertex, z 
per occupied vertex, times N x ^ /\Aut(T) \ for the decorated graph T. 

The calculation of the matrix integral ( |3.6|) goes through a succession of standard 
steps. 

• Step 1 consists in replacing the integral by one over eigenvalues of the matrices, 



using the Itzykson-Zuber formula |16[]. One is left with 



Z(g,z) = J d N ad N bA(a)A(b)e- NT ^ ab -^- 9Z ^- 9 ^ (3.7) 

where a and b are real diagonal matrices of size N x N, and A stands for the Vandermonde 
determinant A (a) = det[af~ }l<i,j<N- 

• Step 2 consists in first noting that A(a) = det[pj_i(aj)]i<j i:? <jv, for any sequence 
of monic polynomials pj of degree j. One then picks the two sets pj and q m to be inserted 
in A (a) and A (b) respectively in such a way that they form a bi-orthogonal set wrt the 
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following bilinear form: 

(f,g) = J dxdye- N ^y-* 2 / 2 -° z * 4 / 4 -°y 4 ^f(x)g(y) (3.8) 
formally defined by using the size-one matrix propagators. Writing 

(Pj,Qm) = 5j,rnhj(g, z) (3.9) 

one then computes 



N-l 



• Step 3 consists in calculating the h/s. One introduces the operators Q± and Q2 of 
multiplication by x and y, expressed respectively on the basis Pj(x) and q m (y). Introducing 
similarly Pi = d/dx and Pi = d/dy, one obtains upon integration by parts the following 
relations 

Q2-Q1- 9*Ql 



p (3-ii) 
-§=Q\-aQl 

where the adjoint is obtained by letting the corresponding operator act on the other poly- 
nomial basis. As a consequence of ( |3.11|) , the Q operators take the form 



QlPn(x) = Pn+l(x) + r n p n _i(x) + S 

nPn — 3 

(3.12) 

<?2<?m(y) = q m +i(y) + (y) + (y) 

Note the presence of only every other term, due to the 7Li symmetry of the potential 
under A — > —A, B — ► —B. The fact that Qi and Qi have a finite range is generic in 
multimatrix models with polynomial potentials. The equations ( |3.11 ) finally turn into 



recursion relations for the various coefficients of ( |3.12| ), which upon elimination allow for 
expressing v n = h n /h n -i, and finally the free energy 

F(g, z) = ^Log Z(g, z) = ± ^(iV - z)Log (3.13) 

• As we will be only interested in the large limit, in which only planar graphs 
are selected (with genus zero free energy fo(g,z) = liniTv^oo F(g, z)), let us perform the 
Step 4, namely the reduction of equations in the large N case. We make the standard 
assumption that all coefficients r n ,h n , etc... tend to smooth functions of x = n/N for 
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n, N — > oo while x remains fixed. This is certainly true in the vicinity of the Gaussian 
model, for small g. We have the limits 

v n -> v(x), r n -> r(x), s n -> (3.14) 

and similarly for the tilded quantities. In the same limit, the Q operators (|3.12| ) become 
simply 

Qi -> a + ra -1 + sa~ 3 

(3.15) 

Q 2 -> r + rr 1 + sr 3 

where the commuting dummy variables a, r are inherited from the corresponding (non- 
commuting) index shift operator acting on the respective polynomial bases. In particular, 
the operator adjoint simply follows from = fa -1 . Plugging ( 3.15 ) into ( |3.11| ) and noting 



that Pi/N = xa~ x + 0(cr~ 2 ), we finally get a set of algebraic equations for the r's and s's 
as well as their tilded counterparts. Defining rescaled versions V = gv, R = gr, S = g 2 s 
and similarly for the tilded quantities, we get 

gx = V -R- 3z(R 2 + S) 
= R-V{l + ZzR) 

= S-zV s (3.16) 
= R - 3VR 

o = s-v 3 



These algebraic equations were given a combinatorial interpretation in [|I7j |L8|], where 
the various functions V, R, S, ... were shown to generate decorated rooted trees. After 
elimination, the system ( |3.16|) reduces to 

W 2 

gx = <p(V) = V(l - 3zV 2 ) - {1 _ 9zV2)2 (3-17) 

For fixed g and z this equation defines upon inversion the unique function V(x) such that 
V(x) = gx + 0(g 2 ) at small g. It encodes the asymptotic properties of the sequence v n , 
hence those of the /i n 's, and finally we have 

fo(g, z) = [ dx(l - x)Log (3.18) 

Jo \ gx J 

where we have identified v(0, z) = x. 
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It is now a simple exercise to find the singularities of the thermodynamic function 
fo(g,z). Noting indeed that 

d V^/o( 9 ,,) = ?WM (3-19) 



dg 2 dg dg \ g 

where V 9jZ is simply the unique solution to g = <p(V giZ ) such that V 9jZ = g+0(g 2 ) (identical 
to V(x = 1)), we see that the singularities of the free energy are those of V 9iZ . Defining 
the critical cosmological constant by 

g c (z) = <p{V c {z)) 

u , , 3 - 20 

where V c (z) is the smallest value of V where a maximum of ip is attained, we immediately 
get the critical behavior fo\ s i ng ~ (<7c(#) ~~ flO 2-7 where 7 = —1/2, as this generically 
corresponds to a square root singularity of V from (g — g c ) ~ (p"(V c )(V — V c ) 2 /2. Eqns. 
( |3.20|) are easily parametrized as 

_ 12^(1 + 3u) 2 _ (l-3u) 4 (l + 10u-15u 2 ) 

Z ~ (l-3u)8 9c[Z) - 12(1 + 3u) 2 {6 - 21) 

This gives a critical curve ending at the point 

25 



8192 



(11^ + 25) < 



64 r- 

g.=g c {z.)= (13V5 -29) (3.22) 
45 

32 

y_=y c (^_) = -(7v / 5-15) 

and where in addition we have (p"(V-) = 0, which implies a susceptibility exponent 7 = 
—1/3, as we now have an equation of the form g — g_ ~ if'"(V-)(V — V-) 3 /6. 

To conclude, we find only one critical value of z at a negative value Z- , corresponding 
to a non-unitary critical point. The value 7 = —1/3 can be carefully traced back to the 
coupling of the Yang-Lee edge sing ularity with c(2, 5) = -22/5 to 2DQGi. This shows 



3 This is in apparent contradiction with the formula ( [2.7| ) which would lead to j s t r = —3/2. 
This is due to the fact that the leading behavior observed in the matrix model solution corresponds 
to another scale set not by the identity operator (leading by integration over the random surfaces 
to the area A = J Idd 2 z or to the conjugated coupling g via the off-critical term (g c — g) J Idd 2 z 
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that no crystallization transition takes place when we couple the hard-particle model to 
ordinary 2DQG. 

Let us however complete the task of computing the free energy. Using the Lagrange 
inversion formula, we may explicitly invert ( pM7| ) and evaluate any function of V 9iZ : 



HVg,z) = J2 9 i<( 
n>l J 



2m (f(w) 



Expressing fo{g,z) = J2 n >i fn(z)g n , and using eqn. ( |3.19| ), we find that f n (z) = 
<p n (z)/((n + l)(n + 2)), where <p n (z) are generated by Log(V 9jZ /g) = Y, n >i l Pn{z)g n . We 
may now use the Lagrange inversion formula ( |3.23| ) for h(V) =Log(V). This gives 

, s 1 f dw 1 ,„ „ 

<Pn(*) = -<t T - 3.24 

n J 2ittw (p{w) n 

This finally leads to the compact formula for the free energy f n (z) for hard particles on 
planar tetravalent fatgraphs with n vertices, and a weight z per particle: 



1 (2n -2l-2 V 



n(n + l)(n + 2) ^— ' V n — 1 

v. 1 / V 1 / l,p>0 x 



X 



Hl+P)<n (3.25) 



X 



2n - I - 2p - 1\ An - 21 - p - 2 
I )\ P 



3 n-l z l+p 



As a side remark, we recover the result for pure tetravalent graphs when z = 0, namely 
f n (z = 0) = 3 n ( 2 ?i n )/(2n(n + l)(n + 2)) = N n /n, N n as in (|2~T2|) (the discrepancy is due to 
the extra rooting in ( |2.12| ) which implies N n = nf n ). More interestingly, if we take z — » oo 



in the gravitational action) but by that with the lowest (negative) conformal dimension, as is the 
case in all non-unitary theories with c(p,q), say with q > p + 1. This operator has conformal 
dimension ho = (1 — (p — q) 2 )/(4pq), and gravitationally dressed dimension Ao = (p — q + l)/(2p) 
from ( |2.8| ). It corresponds to a term (/i c — //) f $?j d 2 2; in the gravitational action, where by a 
dimensional argument we find g c — g ~ (/x c — /x) 1_A °. The exponent 7 observed in the matrix 

model corresponds to the scaling in terms of fi c — /U. Once translated back into the language of the 

2 — 1 

correct cosmological constant g c — g, we have f\sing ~ (a*c — A 1 ) 2 7 ~ {9c — g) 1_A ° = (9c — g) 2 l ' ;i!tr . 
This allows to finally identify 7 = —2/{p + q — 1) for any minimal model, while 7 s t r = (p — g)/p 



follows from (2.7). Here, for 7 = —1/3 we must have p + g = 7, hence p = 2,q = 5. The only 
other possibility would be p = 3 and q = 4, but the corresponding model is unitary and could not 
have a negative critical value of z as in ( 3.22| ). 
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in ( 3.25Q , we should tend to a crystalline groundstate in which vertices of the graphs are 



maximally occupied. We find that only j 2n survive, with the result 



3 n_1 /3n N 
hn = n(2n + l)(2n + 2) { n J (3 ' 26) 



The number (2n)/2n counts the rooted bipartite tetravalent planar graphs |19[ |2(J. An- 
other way of understanding this result, is by looking directly at the matrix integral ( |3.6|) . 
Indeed, performing the change of matrix variables A — > A/z^, B — > Bz^ and redefining 
9^9/ we immediately see that in the limit z — > 00 the matrix integral ( |3.6| ) reduces 
to 

Z(g, 00) = J dAdBe- N ^ AB -^ A *+ B4 » (3.27) 

which clearly enumerates the bipartite tetravalent graphs. 

The lesson to be drawn from this section is that summing over arbitrary tetravalent 
graphs destroys the crystallization transition of the hard-particle model. However, the 
crystalline groundstates which would make this transition possible dominate the sum at 
large z, and involve only bipartite tetravalent graphs. The problem is that these graphs 
have a negligible contribution at all finite values of z. This suggests to reduce the range 
of summation to bipartite graphs from the very beginning, which is the subject of next 
section. 

3.4- Hard objects on bipartite trivalent graphs: matrix model solution 

Simply for pedagogical purposes, we will address here the problem of enumeration 
of hard-particle configurations on trivalent bipartite fat graphs. The case of tetravalent 
graphs was solved in [1]], and presents only technical complications, whereas the final 
results concerning the crystallization transition and its universality class are the same. So 
we choose to concentrate here on the case of hard-triangles coupled to 2DQG. 

We wish to generate bipartite trivalent graphs say with alternating black and white 
vertices, which in addition may be either occupied or empty. To generate the desired 
decorated graphs, we now need four matrices A^A^A^A^ each standing for half-edges 
connected to empty white, occupied black, occupied white and empty black vertices. The 
corresponding matrix integral reads 



Z(g,z) = J dA l dA 2 dA z dA A e- NTvV ^ A ^ A ^ 
V(A U A 2: A 3 , M) = AiA 2 - A 2 A 3 + A 3 A 4 - + 3) - gz (3 + 3] 

19 



(3.28) 



The quadratic form in V(Ai, A 2 , A 3 , A 4 ) has been engineered so as to reproduce the correct 
propagators, namely that only black and white vertices are connected in the Feynman 
diagrams ((AiAj) = if i and j have the same parity) and that two occupied vertices 
exclude one-another ((A 2 A^) — 0). 

We may now repeat the straightforward, though tedious, steps 1-4 of the previous 
section. 

• Step 1 takes us to the eigenvalue integral 

Z(g, z) = J d N a 1 d N a 2 d N a 3 d N a 4 A(a 1 )A(a 4 )e- NTlV{ai > a2 > a:i > a * ) (3.29) 

where the Oj are real diagonal matrices of size N x N. 

• Step 2 goes through as well, by simply replacing the bilinear form ( |3.8|) by 



(f,g) = J dx^dxsdxie-^^'^'^'^fix^gixi) (3.30) 

• Step 3: we define analogously operators Qi,Pi, i = 1,2,3,4, acting on the left 
orthogonal polynomial basis for i = 1,2 and on the right for % = 3,4, and satisfying the 
following system, obtained by integration by parts: 

Pi n rr X - 



(3.31) 



¥ =d Xl V(Q u Q2,Q 3 ,Q4)=Q 2 -gQi 
j± = = d Xa V(Q u Q 2 , Q 3 , Qa) =Qi-Qs- gzQl 
^ = = d Xa V(Q u Q2, Qa, Qa) = Q 4 - Q2 - gzQl 
^ = d XA V(Q u Q 2 , Q 3 , Q 4 ) =Q 3 - gQ\ 

The obvious symmetry of the potential under the interchange X{ <-> x$-i allows to infer 
that Q4 = Q\ and Q3 = Q 2 , so that the main equations reduce simply to 

— = Q2 - gQ\ 

N V2 Wl (3.32) 
Q x = Ql + gzQl 

The corresponding Q-operators are easily expressed on the left basis of orthogonal poly- 
nomials as 

Ql = a + a" V 1 ) + a-M 2 ) + a" V 3 ) 

Q 2 = 0*8W+*- 1 8M+*-* 8 W (3 ' 33) 
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where a is the shift operator acting on the p's as ap n = p n +i and the operators 
are diagonal. The presence of only powers of a spaced by multiples of 3 is a consequence 
of the accidental Z 3 -symmetry of the potential V, namely thet V(ux-\_, 0x2,^x3, 0x4) = 
V(xi, X2, X3, X4) for 00 = . The equations ( |3.32| ) together with the fact that P\/N ~ 
(n/N)a~ 1 +0(a~ 2 ) allow in principle for solving the model for all iV by expressing recursion 
relations for the quantity v n = h n /h n -\. 

• Step 4: We are however only interested here in the planar iV — > 00 limit, in 
which o becomes a commuting dummy variable, and the r^, some functions of x = 
n/N, while n, N —>■ 00. The equations ( |3.32|) translate simply into algebraic equations 
for these functions, as well as for v(x) = limt> n . Upon defining the rescaled quantities 
RW = g^ z 2 r^\ S<® = 8<®/g, = g 2 zs^\ = g b z 2 s^ and V = g 2 zv, we find the 
algebraic system 

= 1 

= - zV 4 

V = S {1) + 2S {1) V (3.34) 
R W = V 2 + z(S^) 2 + 2zS^ 
g 2 z 2 x = zSW-2RW 
By elimination, we are simply left with 

g 2 z 2 x = <p(V) = z {1 J 2V)2 ~ 2V*(1 " 2V2 ) ( 3 - 35 ) 

while the desired planar free energy reads fo(g,z) = f Q (1 — x)Log(V / (g 2 zx)) by virtue 
of the obvious generalization of ( |3.18[ ). The function V must be computed by inverting 
( p.35|) order by order in g 2 x. As in the case of previous section, the singularities of the 
planar free energy come from those of the solution V 9jZ of (|3.35| ) at x = 1, as we still have 
a relation of the form d 2 /dg 2 g 3 d/dgfo(g, z) = gd/dgLog(V g , z /(g 2 z)). 

Expressing the critical points g c (z) solving the obvious adaptation of ( |3.20|) , namely 
that g c (z) 2 = (p(V c (z)) / z 2 , while (p'(V c (z)) = 0, with the appropriate value of (p from 
( |3.35|) , we see a first difference with the case of previous section, namely that as 

<p'(V) ~ (l-2V)(z-4V(l + 2Vf) (3.36) 

we now have two critical curves 

9l{z) = h~^ z>0 ' v = 1/2 

1 1 sy 1 IDT/ 2 (3.37) 

gl{z) = 8(1 +2vr 1 with z[v) = 4V{1 + 2Vf 
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Fig. 6: The critical curves g 2 (z) are represented in solid (g%) and dashed 
(gf) lines. The higher critical endpoint Z- = — 2 9 /5 5 and intersection point 
z + = 32 respectively correspond to the Yang-Lee edge singularity and to the 
critical Ising model both coupled to 2DQG. 



along which we generically have a string susceptibility exponent of 7 = —1/2. These 
are displayed in Fig. 6. While the curve g\ (z) possesses an endpoint at a value z = 
Z- = — 2 9 /5 5 (where V_ = —1/10, and g 2 _ = 3.5 7 /2 20 ) at which (p" vanishes and which 
qualitatively resembles that found in the previous section, we now get another critical 
point at the intersection of the two curves g\(z) and g\{z) at z = z + = 2 5 (with V+ = 1/2 
and g\ = 15/2 12 ) where we also have a vanishing of ip" . Both critical points correspond 
to 7 = —1/3, but while the first one is at a negative value of z (hence expected to be 
described by a non-unitary CFT), the other one is at a positive value of z. The only two 
minimal model candidate CFT's for 7 = — 1/3 are such that p + q = 7 (from the relation 
7 = — 2/(p + q — 1)). This leaves us with only the Yang-Lee edge singularity, non-unitary 
with c(2, 5) = —22/5, and the critical Ising model, unitary with c(3,4) = 1/2. Noting 
that the second critical point found separates two phases of low and high occupancy, we 
deduce from the above discussion that the endpoint corresponds to the Yang-Lee model 
coupled to 2DQG, while the intersection point corresponds to the crystallization transition, 
and is described by the critical Ising model coupled to 2DQG. This assertion is further 
confirmed by actually computing critical exponents corresponding to various operators 
(see ]1J for details). This completes our gravitational proof that the hard triangle model's 
crystallization transition lies in the universality class of the critical Ising model. 

Let us for completeness express the planar free energy's coefficients f2 n {z) of g 2n in 
the (even) g-expansion, by inverting ( |3.35| ) via the Lagrange inversion formula. Repeating 
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the arguments of the previous section, we arrive at 

0<2p<j<n 



w^^t%^ e KJ „-i XYJG-^ (3 ' 38) 



Note this time the presence of signs in the polynomial expression for f2n(z). Note also that 
for z — we recover the number of bipartite trivalent planar graphs with 2n vertices N2 n = 
2nf2n(z = 0) = 2 n ~ 1 ( 2 ™)/((n + l)(n + 2)) dual to the number of Eulerian triangulations 
of area 2n |19| P0|| . The above results also translate into large n asymptotics of f2n{z) via 
( p.l0|) , namely that 



hn{z) ~ rftJ 9c{z) = M Z_<Z<Z + 



g c {z)=gi(z) z>z + 



9c(z ± )- 2n ^ (g e (z-) =v / 15. 



f2n(z±) jo With 



(3.39) 



•2TTT 



^3 l9c{z+) = ^2fi- 

3.5. Partial conclusion 

Using matrix model techniques, we have been able to observe how geometrical con- 
straints (on the hard-particle models) have made bipartiteness of the random lattices 
relevant. Indeed, imposing bipartiteness of the underlying graphs from the very beginning 
proves to be instrumental in restoring the expected crystallization transition. This allows 
to give a formal gravitational proof in the cases of hard triangles (as well as hard squares) 
that this transition is indeed in the universality class of the 2D critical Ising model. Some 
generalizations of the models shown here have been worked out in and concern parti- 
cle models subject to a weaker exclusion constraint. These allow in principle to visit all 
critical points described by minimal models coupled to 2DQG, and must all be defined on 
bipartite graphs as well. 

One could wonder what happens in the case of hard hexagons. Inspired by the previous 
lesson, we would expect that the crystallization transition for hard hexagons only survives 
on tripartite (i.e. vertex-tricolored) fatgraphs. We expect indeed the existence of three 
symmetric competing maximally occupied crystalline groundstates to induce a transition 
in the universality class of the critical 3 state Potts model coupled to 2DQG. Attempts in 
this direction show that a certain rectangular matrix model should do the job of generating 
the corresponding configurations of hard particles on tricolored trivalent graphs, but this 
model has not been amenable to an exact solution yet |22 . 
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As a both final and preliminary remark for the remainder of this note, the crucial 
properties of colorability of the random lattices we have summed over have guaranteed 
the existence of crystalline groundstates, also responsible for the existence of the phase 
transition. In the next example of meanders, we will also arrive at a similar discussion of 
critical properties according to the colorability constraints of the random lattices we will 
sum over. However in this new case, not only the structure of groundstates does rely on 
the colorability constraint, but the very definition of the degrees of freedom of the model 
as well, essential in determining the central charge of the underlying CFT. 



4. Meanders 

In this second part, we address the apparently unrelated problem of enumeration 
of meanders, namely of topologically inequivalent configurations of a road (simple closed 
curve) crossing a river (infinite line) through a given number 2n of bridges (simple inter- 
sections). After first describing the problem in detail, we will turn to a description of the 
meander configurations as some particular random tetravalent planar graphs decorated by 
loops (actually two loops: the road, and the river). This leads us naturally to study the 
more general problem of loop models coupled to 2DQG. In carefully pursuing this, we 
will eventually see that geometrical constraints on these models (the so-called fully-packed 
loop models) make again the colorability of the underlying lattice relevant, like in the 
hard-particle case. We will be able to identify the meander universality class by a subtle 
reasoning on the type of random lattices we must sum over to finally generate the correct 
meander configurations. 



4-1. The meander enumeration problem 

Meanders are defined as planar graph configurations of a closed nonselfintersecting 
curve (road) crossing a line (river) through a fixed number 2n of simple intersection points 
(bridges). These configurations must be counted up to smooth deformations preserving 
the topology. We denote by Mi n the total number of such distinct configurations for given 
number of bridges In. The meander enumeration is an old problem: it can probably be 
traced back to some work by Poincare (1911), and reemerged in various contexts since: as 
mathematical recreation J23|, as folding problem [53| [P5 |, in relation to the 16th Hilbert 
problem [25], in the theory of invariants of 3-manifolds |]2~T ], in computer science J2S 
abstract algebraic terms [p9[ [30], and in its own right |[1] [32] []33| [^4 . 
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In the works [25[] and subsequent, our main motivation was the study of the folding 
problem of polymer chains. Such a polymer is ideally described by a chain of identical 
line segments attached by their ends, which serve as hinges between adjacent segments. 
Think of a single strip of stamps, which can be folded along the edges common to each 
neighboring stamps. We will distinguish between closed and open polymers according to 
whether the chain forms a loop or is open with two free ends. We will be addressing the 
compact self-avoiding folding of such objects, namely study the various ways in which the 
polymer can be completely folded onto one of its segments. Note that a closed compactly 
foldable polymer must have an even number of segments. 



closed 




open 



Fig. 7: Two compactly folded polymers and the corresponding meanders. 
The first example is a closed polymer with 8 segments, and corresponds to 
a meander with 8 bridges. The second example is an open polymer with 6 
segments, and corresponds to a semi-meander with 7 bridges. 

To distinguish between the various ways of compactly folding a closed polymer, we 
will represent the folded objects as a meander with In bridges. To visualize the relation 
between compactly folded closed polymers and meanders, it is simplest to imagine we 
draw a line perpendicular to the segments forming the folded polymer with a total of 
In intersections (each segment intersects the line once), and then separate the various 
segments (see Fig. 7). 

In the case of an open polymer with say n — 1 segments, let us attach one of its ends 
to say a wall or a support (see Fig. 7), so as to prevent the polymer from winding around 
that end (this is exactly the situation in a strip of stamps, attached by one end to its 
support). Starting from a compactly folded configuration, let us again draw this time a 
circle that intersects each of the n segments once, and also intersects the support once. 
Extending the polymer so as to let it form a half-line with origin its free end, we form a 
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planar configuration of a non-selfintersecting loop (road) crossing a half-line (river with 
a source) through n points. These configurations considered up to smooth deformations 
preserving the topology are called semi-meanders (see Fig. 7 for an illustration). The total 
number of semi-meanders with n bridges is denoted by M n . 

The aim of the subsequent sections is to give a semi-rigorous physical argument lead- 
ing to the prediction of the meander and semi-meander configuration exponents a, a 
which govern the large n asymptotics of the corresponding numbers Mm ~ g~ 2n /n a 
and M n ~ g~ n /n a , g c some constant. This is done by identifying the universality class 
of the corresponding critical phenomenon, characteristic of the compact folding of one- 
dimensional objects. We try here to summarize and put in perspective various works 
among which the paper where these results were first obtained. 

4-2. Fixed lattice results for (fully-packed) loop models 



1/K 

A 

+oo 




Fig. 8: Qualitative phase diagram of the O(n) model on the honeycomb 
lattice. The thick lines correspond to critical curves, for which we have in- 
dicated the value of the central charge. We have parametrized n = 2cos7re, 
< e < 1. 
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Loop models have been studied extensively for decades as toy models for describing 
(self-avoiding) polymers and we will only summarize the results obtained for their critical 
behavior in two dimensions. They are defined as statistical lattice models, for which the 
configuration maps a go from the set of edges of the lattice to T = {0, 1}, where cr(e) = 1 
iff the edge e is occupied by a loop edge. 

More specifically let us first discuss the case of loop (so-called O(n)) models on the 
honeycomb lattice. The self-avoidance constraint is very easy to implement as at most 
one loop may visit any given vertex, hence at each vertex with adjacent edges e,e',e", 
we either have cr(e) = cr(e') = cr(e") = 0, or exactly two of the edges are occupied say 
a(e) = cr(e') = 1, while cr(e") = 0. We attach two types of Boltzmann weights to the loop 
configurations: a weight K per occupied edge, and a weight n per loop. The partition 
function of the model reads 

Z(n,K)= J2 nL RE ( 41 ) 

loop configs. 

where L is the total number of loops and E the total number of loop edges. The qualitative 
phase diagram of this model is displayed in Fig. 8. The model undergoes various phase 
transitions within the range —2 < n < 2, in which we parametrize n = 2 cos7re, < e < 1. 
It has three critical curves (represented as a zig-zag thick line and a thick straight line), 
characterized by the density of loops being either 0, between and 1 or 1, respectively 
called dilute, dense and fully-packed phase. The latter case corresponds to a situation 
where the lattice is maximally occupied by loops, namely each vertex is visited by a loop, 
and it is obtained in the limit K — > oo. The three critical curves are respectively described 



by CFT's with respective central charges [[35j |36| 



dilute : c = 1 — 6- 



e 2 



1 + e 



e 2 

dense : c = 1 — 6 (4-2) 

1 — e 



fully — packed : c = 2 — 6 



e 2 



1 



It is beyond the scope of these notes to explain in detail these three formulas. They 
involve the so-called Coulomb gas description of the various critical theories, namely a 
description involving scalar free fields, each contributing 1 to the central charge, while a 
suitably defined background electric charge accounts for the contribution proportional to 
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e 2 . Typically, these scalar fields are the continuum limit of discrete "height variables" 
describing the degrees of freedom of the original model. In the case of loops, one may 
indeed define a dual height variable as follows. We first orient the loops arbitrarily, and 
then define a height variable as constant within each domain delimited by loops, and 
increasing or decreasing by a fixed amount when one goes across a loop pointing to the left 
or right (in other words, the loops form a contour plot of the height variable). This explains 
at least vaguely why only one scalar field will in general be sufficient to give a complete 
description of loop models. The background electric charge is also easy to understand as 
follows. To produce the correct weight n = 2 cos7re per loop, one simply attaches a weight 
e ±me £ Qr eacri orientation. In turn, this is obtained by simply introducing a local weight 
e ±i7re/6 on eacri ver tex visited by a loop, according to whether the loop makes a left or 
right turn. Indeed, on the hexagonal lattice the total of left minus right turns is always 
±6, and the weight n per loop follows from the summation over all loop orientations. Now 
if we try to make sense of the loop model in a cylindric geometry, we see that if a given 
loop winds around the cylinder, it will have as many left and right turns, hence receive a 
wrong weight 2 instead of n. This is repaired by adding the background electric charge, 
which modifies the free field action and restores the correct weight, while modifying the 
central charge accordingly. 

The fully-packed case is clearly distinguished in that it requires two scalar fields, rather 
than just one for the two other cases. More precisely, the central charge for fully-packed 
loops is exactly one more than that of the dense loops for the same value of n. Let us 
show in the two particular cases n = 1, 2 where this extra degree of freedom comes from. 

The case n = 1 describes various equivalent problems for which configurations are 
depicted in Fig. 9: (i) rhombus tiling of the plane (ii) groundstates of the antiferromagnetic 
Ising model on the triangular lattice (iii) fully-packed loops on the honeycomb lattice. The 
interpretation (i) allows to immediately understand why the central charge c = 1 in this 
case. Indeed, the rhombus tiling may also be regarded as a view in perspective of the 
interface between two media, made of the piling- up of cubes. The natural height variable 
for this problem is simply the actual height of the various cubes, which in the continuum 
becomes the z coordinate of some landscape. 

The case n = 2 may be interpreted as the problem of "phantom folding" of the 2D 



triangular lattice ||38|| . We start from the triangular lattice, viewed as a surface, namely 
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Rhombus Tiling 



AF Ising Groundstate 



Fully-Packed Loops 



Fig. 9: Various interpretations of the fully-packed loop model's configura- 
tions on the hexagonal lattice: (i) rhombus tiling, obtained via the dual of the 
hexagonal lattice, the triangular lattice, in which we erase the edges dual to 
unoccupied edges, thus forming rhombi, which come in three possible orien- 
tations, (ii) Antiferromagnetic groundstates of the Ising model: one tries to 
maximize the number of antiferromagnetic interactions (-1 — ) between neigh- 
boring spins on the triangular lattice. Exactly one third of the edges remain 
frustrated (dashed here) namely are adjacent to equal spins: these frustrated 
edges of the triangular lattice are nothing but the erased ones in the rhombus 
tiling, (iii) Fully-packed loops on the hexagonal lattice: each vertex is visited 
by exactly one loop. 





Fig. 10: The triangular lattice is represented as a surface, with edges viewed 
as tangent vectors (orientations are indicated by arrows), such that the sum 
over edges around each face vanishes Xlface ^ = 0- A phantom folding map a 
is simply a map sendig tangent vectors to either of the three images 1,2,3 that 
we have represented. This image is simply the value of the tangent vector in 
the folded configuration. As such, it must satisfy the face-rigidity constraint 
that Eface <KeO = 0. 
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with edges drawn as tangent vectors e, in such a way that the sum of tangent vectors 
around each face vanishes (see Fig. 10), namely 



E 



e = 







(4.3) 



e around face / 



A phantom folding configuration is simply a configuration of this surface after being folded 
back to the plane, in a process where faces are not deformed, and edges serve as hinges 
between adjacent triangles. We are however interested only in the final folded states, not 
in their realizability with an actual material such as a piece of paper. It means that some 
of the foldings we consider may require cutting the sheet and letting it cross itself, hence 
the name phantom. More precisely, a folding map a simply sends each edge e to its new 
value in the folded configuration u{e). Finally the face-rigidity constraint turns into 



It is easy to see that the images a(e) may only take one of three possible values, denoted 
1,2,3 in Fig. 10, namely the three unit vectors with vanishing sum (up to a global ro- 
tation). This allows to replace any phantom folding configuration by an edge tricoloring 
configuration of the triangular lattice, in which we paint edges with colors 1,2,3 in such 
a way that the three edges adjacent to each face have distinct colors. Representing the 
dual of such a configuration, we obtain an edge-tricoloring of the hexagonal lattice. Let us 
concentrate on the edges colored 1 and 2. It is easy to see that edges of alternating colors 
1,2, 1,2, 1,2... form fully-packed loops on the hexagonal lattice. Now exchanging the colors 
1 <-> 2 along any of these loops independently also produces an admissible and distinct 
edge-tricoloring. This shows that the fully-packed loops receive naturally a weight n = 2 
per loop. Finally picking arbitrarily the edges of color 3 in such a way that there be exactly 
one per vertex allows to generate all fully-packed loop configurations. We therefore end up 
with the fully-packed loop model on the hexagonal lattice, with n = 2. The height variable 
of the folding model is easily identified as the actual coordinate h (in IR 2 ) of the vertices 
in the folded configurations. Being two-dimensional, these give naturally rise to two com- 
ponent scalar fields in the continuum, which explains the central charge c = 2 in this case 
(no background charge is necessary as there are no specific left/right turn weights). Note 
that the edge-tricoloring problem of the hexagonal lattice was solved by Baxter [|3!J by 
Bethe Ansatz techniques, and the central charge c = 2 may be directly extracted from 





(4.4) 



e around face / 
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Fig. 11: A typical Fully-Packed loop configuration on the square lattice. 
Assuming doubly periodic boundary conditions, there are 6 black loops (solid 
line) and 4 white ones (dashed lines). Up to rotations, the vertices of the 
model are of the two types (a) "crossing" or (b) "tangent". 



that solution. The general fully-packed result (|4.2| ) was derived by analogous techniques 
in 



37[. 



Let us also mention another interesting loop model, which will eventually allow us 
to solve the meander problem: the two-flavor loop model on the square lattice. The 
configurations of this model are made of loops of two colors (flavors) say 1, 2 (or black, white) 
occupying edges of the square lattice, themselves belonging to at most one loop. The new 
condition is that although these loops are self-avoiding, two loops of distinct colors may 
coexist at a given vertex, by either crossing one-another, or being "tangent to one-another" 
as displayed in Fig. 11 (a) and (b). The model is further defined by attaching a weight 
Ki per edge occupied by a loop of color i = 1,2 and rii per loop of color i = 1, 2. The 
partition function of this model reads: 



Z(n 1 ,n 2 ;K 1 ,K 2 )= ^ n^n^K^K^ (4.5) 

loop configs. 



where Li is the total number of loops of color £ = 1, 2 and Ei the total number of edges 
in the loops of color i = 1,2. In a way similar to the hexagonal lattice case, the model 
undergoes a number of interesting phase transitions, and in particular it still has a dense 
critical line, as well as a fully-packed one, for —2 < ni,ri2 < 2, which we still parametrize by 
rii = 2cos7rei, < e» < 1, £ = 1, 2. The fully-packed constraint restricts the model to only 
configurations where all vertices are visited by the two types of loops, which corresponds 
to taking Ki, K 2 — > oo (see Fig. 11 for an example). This model is referred to as the FPL 2 
model in the following. As before, Coulomb gas techniques have been used to identify 
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precisely the universality classes of the various critical models, leading in particular to 
formulas for the central charge of the underlying CFT EDI : 



dense : c = 2 — 6 



fully — packed : c = 3 — 6 



P 2 


+ 




p 2 


I -ex 
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- e 2 
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1-ei 
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(4.6) 



hence we still have an increase of +1 in the central charge when going from the dense to 
the fully-packed case at fixed m, n^. Again, the 2 and 3 in the central charges correspond 
to the number of degrees of freedom necessitated by the field theoretical description of 
the continuum limit, while the other contributions come from background charges which 
ensure the correct weights rii per loop of color i in cylindric geometries. Let us now show 
where the three scalar fields come from in the fully-packed case. 




-© 
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>- -o C 
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Fig. 12: A typical configuration of the FPL 2 model together with the bi- 
coloration of its vertices (checkerboard of filled (•) and empty (o) dots). We 
have added the corresponding dictionary that allows to map the loop config- 
urations onto A, B, C, D labelings of the edges. 



To do this, let us rephrase the model as a (3D) height model as follows. Starting 
again from an oriented fully-packed black and white loop configuration, we first bicolor 
the vertices of the square lattice, say with alternating filled (•) and empty (o) dots. Then 
we use the dictionary of Fig. 12 to assign one of the four labels A, B, C, D to each colored 
and oriented edge. With this convention, it is clear that edges of type ABAB... alternate 
along black loops, whereas edges of type CDCD... alternate along white loops, and that 
each vertex has one incident edge of each type A, B, C, D. It is clear that the four-labeling 
with A, B, C, D is in one-to-one correspondence with the coloring and orientation of edges 
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Fig. 13: Rules determining the change of the height variable across labeled 
edges. We adopt the Ampere convention that the height is increased (resp. 
decreased) by the edge value if the arrow of the edge points to the left (resp. 
right) . The edge labels must be interpreted as three-dimensional vectors with 
the respective values A, — B, C, — D 

of the FPL model. In particular, the orientation of a given black or white loop is reversed 
if we interchange the A and B or C and D labels along the loop. 

The above colors allow for defining a dual vector height variable on the center of each 
face of the lattice. Indeed, viewing as vectors the A, B, C, D labeling of the edges of the 
lattice, let us arbitrarily fix the height to be zero on a given face of the lattice, and define it 
on all faces by successive use the rules of Fig. 13 for the transition from a face to any of its 
neighbors. Note that it is necessary to impose the condition A + B + C + D = 0to ensure 
that the heights are consistently defined around each vertex. We may therefore assume in 
all generality that A, B, C, D are actually four vectors in M 3 with vanishing sum; let us 
take for definiteness A, B, C, D to be the four unit vectors pointing from the center of a 
tetrahedron towards its vertices. The heights are then clearly three-dimensional, as linear 
combinations of A, B, C, D. 

To conclude this section, we have presented two types of loop models on the hexagonal 
and square lattice and studied their critical lines. In both cases, the fully-packed phase is 
described by one more scalar field than the dense one, resulting in an increase of +1 in 
the CFT's central charge. 

4-3. Fully-packed loop models and 2DQG 

In the next section, we will show how meanders may be viewed as configurations of 
the two-flavor loop model on a random lattice. To pave the road for further results, let us 
first consider here the ordinary (one-flavor) loop model on a random (trivalent) lattice. 

Let us go back for a while to the fully-packed loop model with n = 2, interpreted 
above as the model of phantom folding of the triangular lattice. Recall that in that case 
the height variable has been identified with the position in IR 2 of the vertices in the folded 
configuration. We wish now to define this model on a random triangulation. For this to 
make sense, we need the triangulation to be foldable. Indeed, if it is not, the height variable 
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Non-Eulerian Eulerian 



Fig. 14: Sample configurations of a planar non-Eulerian triangulation, and 
a planar Eulerian one. The Eulerian triangulation has only vertices of even 
valence, which allows to bicolor the triangular faces as shown. 

can no longer be defined as the (2D) position in the folded configuration, but rather like 
in the dense case by viewing the loops as a contour plot for the (ID) height. In genus 
0, the foldability amounts to the fact that each vertex is adjacent to an even number of 
triangles. In other words, the triangulation may be bicolored. Such a triangulation is also 
called Eulerian, as this guarantees the existence of an Eulerian path, visiting all triangles 
of the triangulation. Sample Eulerian and non-Eulerian triangulations are depicted in Fig. 
14 for illustration. Summing over such triangulations will be called Eulerian gravity for 
obvious reasons. Now we may define two possible "gravitational" models of fully-packed 
loops on random trivalent graphs. 

• We may sum over arbitrary trivalent graphs. As the dual triangulation will in 
general not be Eulerian, it will not be foldable, and the extra degree of freedom will be 
lost. The string susceptibility of the corresponding gravitational model (with as usual a 
weight g per trivalent vertex) will be computed using ( |2.7|) with the dense central charge 



e 2 

ordinary gravity : c = 1 — 6 , n = 2cos7re (4.7) 

1 — e 

• We may sum over bipartite trivalent graphs, whose dual triangulation is automat- 
ically foldable, thus preserving the height variable in 1R 2 . The string susceptibility must 
be computed using ( |2.7| ) with the fully-packed central charge 

e 2 

Eulerian gravity : c = 2 — 6- , n = 2 cos ne (4-8) 
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So we expect here a pheomenon very similar to that occurring in the case of hard 
particles. The universality class of the critical point depends crucially on the colorability 
property of the underlying random lattices. This similarity is no coincidence, as the very 
existence of completely folded groundstates relies crucially on the bipartiteness of the lattice 
just like the crystalline groundstates of the hard-particle model did. 

To make the above argument even more solid, let us give a matrix model derivation of 
the results (4~7) and (4T8) in the case n = 1. The fully-packed model on random trivalent 



graphs is described by the partition function 

Z(g) = J dAdBe- NT ^+^- 9AB ^ (4.9) 

where as in Sect. 3.2 we use matrix elements of A and B, two N x N Hermitian matrices, 
to generate empty (resp. occupied) half-edges glued into empty (resp. occupied) edges 
via the propagators (A A) and (BB). Note that the measure of integration in ([O]) is 
normalized in such a way that Z(g = 0) = 1. As now usual, Log Z(g) generates a sum 
over connected trivalent fatgraphs decorated by fully-packed loops of B matrix elements. 
We have n = 1, as there is no extra weight per loop. Noting that the dependence of Z(g) 
on A is Gaussian, we may explicitly integrate over A, with the result: 

Z(g) = [ dBe~ NT ^- 292 ^ (4.10) 



As we already mentioned in Sect. 3.2, this latter integral simply generates random quadran- 
gulations, therefore is a model of pure gravity with no matter, i.e. with c = 0, corresponding 
to the dense formula (|4.7|) , with n = l, e = l/3. 

To generate bipartite trivalent graphs (dual to Eulerian triangulations) , we must use 
complex matrices A, B, whose pictorial representation bears an extra orientation, to distin- 
guish between A and A\ The matrix integral for the fully-packed loop model on bipartite 



trivalent graphs reads [41 



Z E (g) = J dAdB e- NT ^ AAt + BBt -^ AB2 + A ^ Bt ^ (4.11) 

where A and B are complex N x N matrices, and the measure is normalized in such a 
way that Z E {g = 0) = 1. We note that the dependence on A is still Gaussian, hence upon 
integrating it out, we get 

Z E (g) = J dBe - N ^ BB '-^ B2{ - B '^ (4.12) 
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This turns out to be a particular case of the six-vertex model coupled to 2DQG solved 
in 



[Eq| . The corresponding central charge is c = 1, thus confirming the formula ( f4.8|) at 
n = 1, e = 1/3. 

The lesson to be drawn from this section is simple: the extra degree of freedom of the 
fully-packed loop model is wiped out if we are not careful enough with the type of random 
lattices we sum over. Only if these are bipartite will the height variable survive and 
lead eventually to the central charge ( f4.8|) . Otherwise, the mismatch between the random 
lattices and the model takes the universality class back to that of the dense model. We 
expect this result to be quite general and we will apply it in the next section to the case 
of the two-flavor fully-packed loop model on random tetravalent graphs. 

4-4- From fully-packed loops on the square lattice to meanders 

We are now ready to turn to the meander problem. It is by now clear that a meander 
configuration is nothing but a configuration of the two-flavor loop model on a random 
tetravalent graph, with in addition exactly one loop of each color, say the river (color 1) 
and the road (color 2). This assumes that we have closed the river (infinite line) into a loop 
as well. The operation of cutting the river loop and extending it into a line is equivalent 
to a marking (rooting) of the corresponding diagram, and yields an overall factor of 4n 
for diagrams with 2n vertices. In the language of the two-flavor fully-packed loop model, 
extracting the configurations with one loop of each color corresponds to letting ni, n% — > 0, 
in order to select the leading contribution to the partition function, proportional to n±U2 
in this limit. But why at all consider the FPL 2 model? Indeed, only the "crossing" vertex 
of Fig. 11 (a) occurs in a meander configuration. The problem is that a two-loop flavor 
model on the square lattice with only the crossing vertex allowed is trivial, as loops should 
all be straight lines (closed into loops say by periodic boundary conditions). It simply 
means that the degrees of freedom of the lattice version of the meander model are not 
taken into account by this naive definition. We will show that the FPL 2 model provides 
the correct lattice version of the meander model. 

Let us now study the FPL 2 model coupled to 2DQG, namely defined on random 
tetravalent graphs. As explained in Sect. 2, the genus zero partition function for this 
model reads 

Z(n 1 ,n 2 ;g,x,y)= £ 1 J] n^n^(gx) v ^ v \gy) v ^ (4.13) 

tetravalent planar ' V / I FPL configs. 

graphs r on T 
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Fig. 15: A sample bipartite tangent meander with 8 vertices, among which 
6 are crossing vertices (bridges) and 2 are tangency vertices. The vertices are 
bicolored, and all river or road edges connect only vertices of distinct colors. 



where the sum extends over all planar four-valent graphs T, and |Aut(T)| is the order 
of the symmetry group of T, while we have also denoted by V a , V& the total numbers of 
vertices of type a and b defined in Fig. 11 in the particular loop configuration, namely we 
have weighted each crossing of a black and a white loop by gx and each tangency by gy. 
From the lesson of Sect. 4.3, we know that if we wish to preserve all degrees of freedom of 
the model, we must couple it to Eulerian gravity, namely consider it on bipartite random 
tetravalent graphs in the sum (|4.13|) (we will distinguish this particular summation by an 
index E for Eulerian in the notation for the partition function). Indeed, recall that the 
definition of the (3D) height variable in Sect. 4.2 has relied crucially on the preliminary 
bicoloration of vertices of the lattice, allowing for defining the edge variables A, B : C, D. We 
may therefore use ( |4. 13| ) to generate "bipartite tangent meanders" , namely configurations 
of a closed road crossing or tangent to a river (infinite line) through a given number of 
bicolored points 2n, such that only points of distinct color are connected by river or road 
edges. Let \iin denote the number of bipartite tangent meanders with 2n vertices (see Fig. 
15 for an illustration). These are well enumerated by counting configurations of the FPL 2 
model at n\ = n-i = 0, namely e\ = e<i = 1/2, and by taking x = y = 1 in ( 4.13 ). More 
precisely, we have 

hm J-(Z E (m, n 2 - g, 1, 1) - 1) = £ (4.14) 

ni,n 2 — >0 71177,2 

n>l 

where we have expressed that the opening of the river into an infinite line amounts to the 
marking of a river edge, and results in a factor 4n between the fixed n free energy f^ n an d 
the tangent meander number \i% n . Plugging the values e\ = e% = 1/2 into the formula 
(fO|) and using (pT7j), we find 

c=-3 7str = -l±^I (4.15) 
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The free energy at fixed area A for this model is nothing but 
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by directly applying ( |2.10| ). Note that the value of g c is not predicted by this argument, 
only the configuration exponent. We finally deduce the large n asymptotics of the bipartite 
tangent meander numbers 
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(a) (b) 

Fig. 16: Typical height configurations around (a) a crossing vertex and (b) 
a tangency vertex. When going diagonally say from the NW to the SE face 
(shaded), the height increases by an amount restricted to ±(A — C), ±(B + C) 
in the case (a) and only to ±(A + B) in the case (b). Forbiding (b) will 
therefore amount to a reduction of the height variable range from 3 to 2 
dimensions. 

This is not however the end of the story. Indeed, the numbers \i2 n might have noth- 
ing to do with the meander numbers M<i n . We must examine more closely the role of 
the tangency points. Note that the meander numbers are correctly generated by taking 
ni,n,2,y -> in Q4.13|) : 

lim J_ (ZE(ni5 n2 . g , 1, 0) - 1) = V (4.18) 

ni,n 2 — >0 77,177,2 4n 

n>l 

which just amounts to forbiding the tangency vertex b of Fig. 11. Forbiding the tangency 
vertex b will further restrict the range of the height variable of the model. In Fig. 16, 
we have displayed typical height configurations around a type a and a type b vertex for 
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comparison. It is easy to see that the increase in height when going diagonally from a 
face to another through the vertex (shaded faces in Fig. 16), there is more freedom at the 
crossing vertex a, where we may have h — > h± (A — C) or h ± (B + C) , while at the tangent 
vertex b we may only go in one spatial direction h — > h ± (A + B). Forbiding the b vertex 
therefore has the net effect of dimensionally reducing the range of the height variable to 
only two dimensions rather than three in the full model. Forbiding tangencies therefore 
takes the central charge back to the dense value, which for e\ = e% = 1/2 reads 



7 „ ir = _£±^H (4. 19) 



Finally, using Q4.18| ) and applying ( |2.10|) , we get the following meander number asymptotics 
for large n 

M 2 „~£^ a= 29Wl46 

n a 12 K J 

To make the above argument more solid, let us present an alternative route lead- 
ing to the result ( f4.19|) . We may have considered in ( }4.13| ) a summation over arbitrary 



(non-necessarily bipartite) tetravalent graphs. In this case, the bicoloration of vertices 
disappears, and we have therefore no way of distinguishing the vector labels A and — B on 
one hand and C and — D on the other hand. This amounts to imposing an extra condition 

A + B = C + D = (4.21) 

This is precisely the condition one would obtain by considering the dense two-flavor loop 
model on the square lattice. Indeed, we still have the dictionary of Fig. 12 for edge labels, 
but as more vertices are allowed, we get some extra conditions for the height variable to be 
well-defined. If we consider a vertex crossed by only a black or white loop we immediately 
see that ( |4.21| ) must be satisfied. This explains the reduction 3 — > 2 in the expression for 



the dense central charge (fO|) compared to the fully-packed one, as the (2D) height is now 
a linear combination of A, C which we may choose to form an orthonormal basis of 1R 2 . 
Hence the simple fact that we sum over arbitrary tetravalent graphs takes the universality 
class of the FPL 2 model back to that of the dense one. For e± = e<i = 1/2, we find precisely 
the result (fyj). 

We must now see where meanders fit into this picture. Recall that we still have at this 
point crossing and tangent vertices, but we have relaxed the constraint of bicolorability. 
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Fig. 17: A type b vertex of the FPL 2 model coupled to ordinary 2DQG, 
together with its dual height configuration. We note that the NE and SW 
heights are identical. We may therefore undo the vertex as shown, which 
explains its irrelevance. 

We must again examine more closely the role of the tangency points. Forbiding them will 
yield the meander configurations as 

lim —(Z(ni, n 2 - g, 1,0)-1) = ]T ^V" (4.22) 

n>l 

Examining the height variable of the dense model around a tangency vertex, we see that 
the height variable is the same in the faces lying outside of the two loops (see Fig. 17). 
This means simply that the tangency vertex is irrelevant in the ordinary gravity situation, 
namely that its presence or absence, although it might affect the precise value of g c , doesn't 
affect the universality class, hence the configuration exponent is independent of y at x = 1 
in ( f4.13| ). We therefore end up with the meander universality class, determined by Q4.21|) . 



Some important numerical checks of the results of this section have been carried out 



by use of Jensen's transfer matrix method [43], which consists in generating the meander or 



tangent meander or bipartite tangent meander configurations by applying a local "transfer 
matrix" along the river, which implements the addition of a new vertex to a previous 
state. The method is powerful enough to give very accurate values for the extrapolated 
configuration exponents, which all confirm the various results obtained above |]44| . Note 
that these are somewhat contested in ||45|| . however the discrepancy with the predicted 
exponents is extremely small and possibly due to subtleties with the large n corrections. 

4-5. The meander universality class: asymptotics of meandric numbers 

We have now identified the CFT underlying meanders, as the dense two-flavor loop 
model with n\ = n 2 = coupled to ordinary 2DQG. The complete knowledge of the 
conformal operator content of this CFT via the Coulomb gas picture gives access to a host 
of meandric numbers which we describe now. 
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The important operators for our present purpose are those identified as generating 
oriented river vertices, namely the operators <pk (resp. <p-k), k = 1, 2, ... which correspond 
to the insertion of a /c-valent source (resp. sink) vertex from (resp. to) which k oriented 
river edges originate (resp. terminate). Using these, we may generate configurations in 
which the river may itself form a complicated though connected oriented planar graph 
with possibly loops and endpoints, either sources or sinks according to the operators. In 
the Coulomb gas picture, these operators create "magnetic" defect lines along which the 



height variable has discontinuities. The operator 0^ has conformal dimension [[10 



k 2 -4 

hk = -gg- (4.23) 

k = ±1, ±2, ... When coupled to 2DQG, these operators get dressed (into ^k) and acquire 
the dimension ( |2.8| ): 

W8 + 3k 2 - Vb 

A k = 2 p= 4.24 

V29-V5 



As a preliminary remark, we note that h±2 = A±2 = 0. ^±2 indeed corresponds 
to the marking of an edge of the river in meanders, and moreover such operators must 
go by source/sink pairs for the orientations of the pieces of river connecting them to be 



compatible. Applying (|2.11|) to the two-point correlator (^2^-2) A at fixed large area A, 
we find 

(^-2) A - (4.25) 

while the meander counterpart (with a closed river) behaves as M A /(2A) ~ g- A /(A 3 ~^). 
We see that the net effect of the insertion of the operators ^±2 is an overall factor pro- 
portional to A 2 , which confirms their interpretation as marking operators. These are very 
important to keep in mind, as they might be required to ensure source/sink balance in 
various river geometries. 

We may now turn to the case of semi-meanders, namely meanders in which the river 
is a semi-infinite line around the origin of which the road may freely wind. Considering the 
point at infinity on the river as just another point, the semi-meanders may equivalently 
be viewed as meanders whose river is made of a segment. Sending one of the ends of 
the segment to infinity just resolves the winding ambiguities around both ends. Using 
the above river insertion operators, we immediately identify the generating function for 
semi-meanders as 

(*!*_!> = J2 M n9 n (4-26) 

n>l 
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Using again ( 2.11 ) and the explicit values of A±i via ( |4.24j ), we arrive at the large n 
asymptotics 



M n ~ 



9c 



a 



I + 2A, %„ - 1 + ^(v^ + v^9) 



(4.27) 



Note that we expect the value of g c to be the same for meanders and semi-meanders, as 
both objects occur as thermodynamic quantities in the same effective field theory. 




Star 



Eight 



Cherry 



Fig. 18: Three types of meandric configurations in which the river has the 
geometry of (a) a /c-valent star (b) an "eight" (c) a "cherry". The vertices 
corresponding to river sinks or sources are represented by filled circles (•). 
The edges of river inbetween them are oriented accordingly. The road (dashed 
line) may freely wind around univalent vertices. 



We may now generate many more meandric numbers by considering more general 
correlators. To name a few (all depicted in Fig. 18), we may generate rivers with the 
geometry of a star with one /c-valent source vertex and k univalent sink vertices generated 
by (\l/fc(\E , _i) A: ), rivers with the geometry of an "eight" with one tetravalent source vertex 
and two loops, each containing a bivalent sink vertex generated by (\&4 (\E r _2) 2 ) , or rivers 
with the geometry of a "cherry" with one trivalent source vertex, one univalent sink, and 
one loop, marked by a bivalent sink vertex generated by (^s^ -2), etc ... For each of 
these situations, we get the corresponding configuration exponent a = 3 — 7 + ^(Aj — 1) 
by applying ( |2 . 1 1| ) with the dimensions Q4.24Q . We get respectively 



«k-star 



height 



-^(v^ + v / 29)(v / 3A^+8 + k(Vll - 2V29) + 4V29 - 2^5) 



48 



i-(V5 + V29)(Vl4 + V5) 



24 



(4.28) 



cherry = ^(v^ + V29)(VU + ^35) 



These values were also checked numerically in 44 . 
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5. Conclusion 

In these notes, we have tried to clarify the role of random lattices when coupling 
geometrically constrained systems such as hard particles or fully-packed loops to 2DQG. 
We have in both cases shown that the application of the famous KPZ formulas ( |2.7|)(pTBD 
can be subtle, and involves the correct understanding of the models' degrees of freedom 
both on the fixed and random lattices. The main lesson is that if we want to preserve the 
essence (universality class) of a given model when coupled to 2DQG, then we must pay 
the price of restricting the fluctuations of space (random lattices) to those and only those 
preserving these degrees of freedom. 

Provided we incorporate this lesson, we are free to use back and forth the KPZ formulas 
to transpose fixed lattice results into random lattice ones and vice versa. In the case of 
hard objects, we have used the 2DQG picture to infer some result on fixed lattices. In 
the case of fully-packed loops, we have used the opposite strategy to solve the meander 
asymptotics problem by using fixed lattice results on the dense or fully-packed two-flavor 
loop model. Along the way we have presented a number of other relations between fixed 
and random lattice models. 

One could still wonder, now that we know the exact answers, whether one could devise 
some more rigorous mathematical proof of the various results inferred. In particular, the 
hard-triangle and hard-square models remain to be solved. As to meanders, the existence 
of a transfer matrix formulation may give some hope to be able to attain a combinatorial 
proof of our assertions. One could also wonder whether the recent combinatorial techniques 
for tackling planar graphs by use of trees might apply to meandric or loop problems (see 
g(J H 0] @ for instance). 



-• — o — • — 6 — • — o — * — 6 — 4 — o- 



Fig. 19: A configuration of the one-flavor fully-packed loop model on a 
bipartite trivalent planar graph. We have opened the loop into a straight 
line and represented the unoccupied edges as dashed non-intersecting arches. 
Both loop and non-loop edges connect only vertices of distinct colors. 
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As a final question, let us mention the following simpler enumeration problem for 
which we also know the exact configuration exponent. It corresponds to the one-flavor 
fully-packed loop model coupled to Eulerian 2DQG, when n — > [^7f. The corresponding 
configurations are simply trivalent planar graphs with bicolored vertices and exactly one 
loop visiting all vertices. Upon opening the loop into an infinite line, we must enumerate 
planar configurations of non-intersecting arches connecting bicolored vertices on the line 
by pairs, in such a way that the colors of the connected vertices are always distinct (both 
along the river and along arches). A sample such configuration is depicted in Fig. 19. The 
resulting number of configurations for a fixed number 2n of vertices, i/2 n , counts also by 
duality the number of Eulerian triangulations with a Hamiltonian cycle. This number has 
the following large n asymptotics 



9c zn 5 _ 13 + V13 
n s ' 6 



(5.1) 



for some constant g c . This is obtained by applying Q2.10 ) with the string susceptibility 
( j2 . 7| ) of the fully-packed loop model with n = coupled to Eulerian gravity, namely with 
the fully-packed central charge ( |4T2] ) with e = 1/2, hence c = — 1, and 7 s t r = — 
Challenge: prove eqn. (pTT|) ! 
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